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Abstract 

We present the results of our QCD analysis for nonsinglet unpolarized quark distributions and 
structure function F2{x,Q'^). New parameterizations are derived for the nonsinglet quark distribu- 
tions for the kinematic wide range of x and Q^. The analysis is based on the Jacobi polynomials 
expansion of the structure function. The higher twist contributions of proton and deuteron struc- 
ture function are obtained in the large x region. Our calculations for nonsinglet unpolarized quark 
distribution functions based on the Jacobi polynomials method are in good agreement with the 
other theoretical models. The values of Aqcd and as{M^) are determined. 
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I. INTRODUCTION 



The deep-inelastic lepton-nucleon scattering is the source of important information about 
the nucleons structure. New and very precise data on nucleon structure functions have 
had a profound impact on our knowledge of parton distributions, in the small and large x 
region. During the last years the accuracy of the obtained experimental data has extensively 
grown up enough to study in detail the status of the comparison of the available data with 
the theoretical predictions of quantum chromodynamics (QCD) in the different regions of 
momentum transfer. 

The importance of deep-inelastic scattering (DIS) for QCD goes well beyond the measure- 
ment of {li- III the past it played a crucial role in establishing the reality of quarks and 
gluons as partons and in promoting QCD as the theory of strong interactions. Nowadays it 
still generates challenges to QCD as, for example, in the domain of structure functions at 
small X 0, 3] or of polarized structure functions or of generalized parton densities jsj and 
so on. 

All calculations of high energy processes with initial hadrons, whether within the stan- 
dard model or exploring new physics, require parton distribution functions (PDF's) as an 
essential input. The reliability of these calculations, which underpins both future theoretical 
and experimental progress, depends on understanding the uncertainties of the PDF's. The 
assessment of PDF's, their uncertainties and extrapolation to the kinematics relevant for 
future colliders such as the LHC is an important challenge to high energy physics in recent 
years. 

The PDF's are derived from global analysis of experimental data from a wide range of 
hard processes in the framework of perturbative QCD. In this work this important problem 
is studied with the help of the method of the structure function reconstruction over their 
Mellin moments, which is based on the expansion of the structure function in terms of Jacobi 



pol yno mials. 
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his method was developed and applied for QCD analysis 
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0, 0, y, y, 0, n, 



16]. The same method has also been applied in polarized case in Refs. 17] 



2l|. 



In this paper we use the deep-inelastic world data for nonsinglet QCD analysis to obtain 
the parton distribution function up to next-to-next-to-leading order (NNLO) approxima- 
tions. The results of the present analysis is based on the Jacobi polynomials expansion of 



2 



the nonsinglet structure function. 

The plan of the paper is to give an introduction of the Jacobi polynomials approach in 
Sec. II. The method of the QCD analysis of nonsinglet structure function, based on Jacobi 
polynomials are written down in this section. In Section III we present a brief review of the 
theoretical formalism of the QCD analysis. A description of the procedure of the QCD fit 
of F2 data are illustrated in Sec. IV. Section V contains final results of the QCD analysis. 
Our conclusions are summarized in Sec. VI. 



II. JACOBI POLYNOMIALS APPROACH 

The evolution equations allow one to calculate the Q^-dependence of the parton distri- 
butions provided at a certain reference point Q^. These distributions are usually parame- 
terized on the basis of plausible theoretical assumptions concerning their behavior near the 
end points x = 0, 1. 

One of the simplest and fastest possibilities in the structure function reconstruction from 
the QCD predictions for its Mellin moments is Jacobi polynomials expansion. The Jacobi 
polynomials are especially suitable for this purpose since they allow one to factor out an 
essential part of the x-dependence of the structure function into the weight function j^. 
Thus, given the Jacobi moments a„,(Q^), a structure function /(x, Q^) may be reconstructed 



in a form of the series 
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max 



Xf{x, Q') = X^{1 - xT <^n{Q^W/{^). (1) 



n=0 



where N^ax is the number of polynomials and G°'^(x) are the Jacobi polynomials of order 

n 

0n'^(^) = E4"^("'/^)^'' (2) 

i=o 

where cj"^(a,/5) are the coefficients expressed through F— functions and satisfy the orthog- 
onality relation with the weight — x)" as in the following: 



dx x^(l - x)"e^'''(x)er^(x) = 5u,i . (3) 
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For the moment, we note that the dependence is entirely contained in the Jacobi moments 

Jo 

n „i 



V / dxx'+\f\a,(3)f{x,Q^' 



= 5^cf)(«,/3)/(j + 2,Q2), 

obtained by inverting Eq. ([1]), using Eqs. ([21 [3]) and also definition of moments, /{j^Q"^) = 
jldx x^-'fix,Q'). 

Using Eqs. ([T]|4]) now, one can relate the structure function with its Mellin moments 

^ max 

n=Q 

n 

xJ2cf\a,mij + 2,Q'), (5) 

j=0 

where F2{j + 2,Q^) are the moments determined in the next section. Nmaxi ol and /3 have 
to be chosen so as to achieve the fastest convergence of the series on the right-hand side of 
Eq. ([5]) and to reconstruct xg\ with the required accuracy. In our analysis we use N^ax = 9, 
a = 3.0 and (3 = 0.5. The same method has been applied to calculate the nonsinglet structure 



unction xF^ from their moments 1^, [l^, ll5| and for polarized structure function xgi 

Obviously the Q^-dependence of the polarized structure function is defined by the 
dependence of the moments. 



III. THEORETICAL FORMALISM OF THE QCD ANALYSIS 



In the common MS factorization scheme the relevant F-?, structure function as extracted 



from the DIS ep process can be, up to NNLO, written as 22, 
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(6) 
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The nonsinglet structure function F2^ns{x,Q'^) for three active (hght) flavors has the repre- 
sentation 



1 



X 



^2,Ns(a;,Q') = C2,NS{X,Q^) ® 



,2^(2)+ 



ix,Q' 



(7) 



The flavor singlet and gluon contributions in Eq. ([6]) reads 



X 9 



.(1) , 0^(2) 



(8) 



^i^2,g(a;,Q') = '^C2,g®g{x,Q^) 



^^2,g + '^^^2,3 



g{x,Q') . 



(9) 



The symbol ® denotes the Mellin convolution 

[A^B]{x)= / (ia;i / (is2 5(x — 0:1X2) A(a;i)i?(x2) 
Jo Jo 



(10) 



In Eq. (CD = m+m — ((i+J) = u^ — dy and = -u+M+rf+J— 2(s + s) = ■u„+(it,+ 2^+2(1— 4s, 
where s = s. Also in Eq. ([8]) S(x, Q^) = Sg=u_rf_s(g + q) = + + 2u + 2d + 2s. Notice 
that in the above equations a = a{Q^) = as{Q^)/4:H denotes the strong coupling constant 
and Cij{N) are the Wilson coefficients [261]. 

The combinations of parton densities in the nonsinglet regime and the valence region 



X > 0.3 for F| in LO is 



-Fiix,Q') 



X 



(x,g^) + -s(x,g^) 



(11) 



where since sea quarks can be neglected 

in the region x > 0.3. So in the x-space we have 

F|(x,Q2)= (^^a;g^g8 + ^xg^g3^ {x,Q^) = ^xu^{x,Q^) + ^x dy{x,Q^) . (12) 

In the above region the combinations of parton densities for are also given by 

F^{x,Q') = (^-^xq+s,s^ (x,Q2) = ^^^u, + d,)ix,Q^) , (13) 



where d = {p + n)/2 and q^^ 3 

In the region x < 0.3 for the difference of the proton and deuteron data we use 



where now g^s : 



2{F^-F,'){x,Q') 

1 1 2 - 

3 ^?NS,3(^'<5^) = 2 - + -x{u- d){x,Q'^) , (14) 

- + 2{u — d) since sea quarks cannot be neglected for x smaller 



than about 0.3. In our calculation we supposed the d — u distribution 

x{d-u){x, Ql) = 1.195a;^-2^(l - xf-^\l + U.05x - 45.52x^) 



(15) 



at Qq = 4 GeV^ which gives a good description of the Drell-Yan dimuon production data 



271]. In our analysis we used the above distribution for considering the symmetry breaking 



of sea quarks [28|, |29|]. By using the solution of the nonsinglet evolution equation for the 
parton densities to 3— loop order [scl], the nonsinglet structure functions are given by 

F,\N, Q') = (l + a C^]Un) + a' C^'Un)) x F,'{N, QI) i^- 



1 ~ ~ '^0) 

Po 



(a^ - al) 



2Po 



PiHN) - ^Po(iV) 
Po 



Po 



^[wrfJ '''''' 



+ 2^ (« - «o)' (A^(iV) - |Po(iV) 



(16) 



Here k = p,d and NS denotes the three above cases, i.e. proton, deuteron and nonsinglet 
structure function. c!f^g{N) are the nonsinglet Wilson coefficients in 0(a™) which can be 
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31 



32| and Pm denote also the Mellin transforms of the (m + 1)— loop splitting 



found in 
functions. 

The strong coupling constant plays a more central role in the present paper to the 
evolution of parton densities. At N'^LO the scale dependence of is given by 

da. 



dlnQ2 



(17) 



fe=0 



The expansion coefficients fik of the /3-function of QCD are known up to A; = 2, i.e., N^LO 



33, 



3^ 



f5o = 11 - 2/3 Uf, 

pi = 102 - 38/3 Uf , 

p2 = 2857/2 - 5033/18 nj + 325/54 



n 



f ' 



(18) 



here Uf stands for the number of effectively massless quark flavors. The strong coupling 
constant up to NNLO is as foUowings |35|: 



1 



1 



jbi ItiLa + 
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(PoLa) 



PqLa 

where La = ln{Q'^ /h?), bk = Pk/Po, and A is the QCD scale parameter. 



[hi [Iu'^La - InL^ - 1) + 62] , (19) 



IV. THE PROCEDURE OF THE QCD FITS OF F2 DATA 

In the present analysis we choose the following parametrization for the valence quark 
densities 

xUy{x,Ql) =Mu a;''"(l - x)''"(l + c„a/x + rf„ x) , 
xd,{x,Ql)=Mdx''^{l-xY'^{l + Cd^ + dd x) , (20) 

in the input scale of Qo = 4 GeV^ and the normalizations Mu and Md being fixed by 
Uydx = 2 and jj^ d^dx = 1, respectively. By QCD fits of the world data for -Ff ''^j we can 
extract valence quark densities using the Jacobi polynomials method. For the nonsinglet 
QCD analysis presented in this paper we use the structure function data measured in charged 
lepton proton and deuteron deep- inelastic scattering. T le experiments contributing to the 
statistics are BCDMS [36], SLAC StI, NMC 38], HI 39|, and ZEUS fl. In our QCD 
analysis we use three data samples : FKx,^^), F^(a;, Q^) in the nonsinglet regime and the 
valence quark region x > 0.3 and -^2^'^ = 2{F2 — -F^) in the region x < 0.3. 

The valence quark region may be parameterized by the nonsinglet combinations of parton 
distributions, which are expressed through the parton distributions of valence quarks. Only 
data with > A GeV^ were included in the analysis and a cut in the hadronic mass of 
W"^ = (i — 1) + rn^ > 12.5 GeV^ was applied in order to widely eliminate higher twist 
(HT) effects from the data samples. After these cuts we are left with 762 data points, 322 
for 232 for F^, and 208 for F^^ . By considering the additional cuts on the BCDMS 



{y > 0.35) and on the NMC data((5^ > 8 GeV^) the total number of data points available 
for the analysis reduce from 762 to 551. 

The simplest possible choice for the function would be 

^data / pi^^^jj rptheor\2 

where AF2f"' is the error associated with data point i. Through ^2^'^°'", is a function 
of the theory parameters. Minimization of would identify parameter values for which 
the theory fits the data. However, the simple form is appropriate only for the ideal case 
of a uniform data set with uncorrelated errors. For data used in the global analysis, most 
experiments combine various systematic errors into one effective error for each data point, 
along with the statistical error. Then, in addition, the fully correlated normalization error 
of the experiment is usually specified separately. For this reason, it is natural to adopt the 



following definition for the effective 4l| : 



Xglobal 

= X1^"^"' labels the different experiments) (22) 



2 



An 



2A 



For the n*^ experiment, AF^^*", and F^^""^ denote the data value, measurement 

uncertainty (statistical and systematic combined), and theoretical value for the i*^ data 
point. AA/'n is the experimental normalization uncertainty and Mn is an overall normalization 
factor for the data of experiment n. The factor Wn is a possible weighting factor(with 
default value— 1). However, we allowed for a relative normalization shift Mn between the 
different data sets within the normalization uncertainties AA/'n quoted by the experiments. 
For example the normalization uncertainty of the NMC (combined) data is estimated to be 
2.5%. The normalization shifts Mn were fitted once and then kept fixed. 

The number of data points for the nonsinglet QCD analysis with their x and ranges, 
and the normalization shifts determined are summarized in Table. I. In this table the first 
column gives (in parentheses) the beam momentum in GeV of the respective data set (num- 
ber), a flag whether the data come from a combined analysis of all beam momenta (comb) or 
whether the data are taken at high momentum transfer (hQ2). The x and range indicate 
in the second and third columns, respectively. The fourth column (F2) contains the number 
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of data points according to the cuts: > 4 GeV^, W"^ > 12.5 GeV^, x > 0.3 for and 
F2 and X < 0.3 for ^2^"^. The reduction of the number of data points by the additional cuts 
on the BCDMS data {y > 0.3) and on the NMC data (Q^ > 8 GeV^) are given in the fifth 
column {F2 cuts). The last column (A/") contains the normalization shifts. 

Now the sums in Xgiobai '^^^^ data sets and in each data set over all data points. The 
minimization of the above value to determine the best parametrization of the unpolarized 



parton distributions is done using the program MINUIT 42|. 



The one a error for the parton density fg as given by Gaussian error propagation is 

-(/.w)^-EE(t)(i)-fe.p.). ^''^ 

where the sum runs over all fitted parameters. The functions dfg/dpi are the derivatives of 
fq with respect to the fit parameter pi, and cov{pi,pj) are the elements of the covariance 
matrix. The derivatives dfg/dpi can be calculated analytically at the input scale Q^. Their 
values at are given by evolution which is performed in Mellin-A^ space. 



V. RESULTS 

In the QCD analysis of the present paper we used three data sets: the structure functions 
Q^) and Q"^) in the region of x > 0.3 and the combination of these structure 

functions F|^^(x,(5^) in the region of a; < 0.3 . Notice that we take into account the cuts 

> 4 GeV^, W"^ > 12.5 GeV^ for our QCD fits to determine some unknown parameters. 
In Fig.([l]) the proton data for F2{x, Q^) are shown in the nonsinglet regime and the valence 
quark region x > 0.3 indicating the above cuts by a vertical dashed line. The solid lines 
correspond to the NNLO QCD fit. 

Now, it is possible to take into account the target mass effects in our calculations. The 
perturbative form of the moments is derived under the assumption that the mass of the 
target hadron is zero (in the limit 00). At intermediate and low this assumption 

will begin to break down and the moments will be subject to potentially significant power 
corrections, of order O {mj^/Q'^), where is the mass of the nucleon. These are known 
as target mass corrections (TMCs) and when included, the moments of flavor nonsinglet 



structure function have the form 



43 



44| 
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Experiment 


X 




GeV^ 




F2 cuts 


A/" 


BCDMS (100) 


0.35 - 0.75 


11.75 


- 75.00 


51 


29 


1.005 


BCDMS (120) 


0.35 - 0.75 


13.25 


- 75.00 


59 


32 


0.998 


BCDMS (200) 


0.35 - 0.75 


32.50 


- 137.50 


50 


28 


0.998 


BCDMS (280) 


0.35 - 0.75 


43.00 


- 230.00 


49 


26 


0.998 


NMC (comb) 


0.35 - 0.50 


7.00 


- 65.00 


15 


14 


1.000 


SLAC (comb) 


0.30 - 0.62 


7.30 


- 21.39 


57 


57 


1.013 


HI (hQ2) 


0.40 - 0.65 


200 - 


- 30000 


26 


26 


1.020 


ZEUS (hQ2) 


0.40 - 0.65 


650 - 


- 30000 


15 


15 


1.007 


proton 






322 


227 





(a) Number of F2 data points. 



Experiment 


X 


GeV^ 


Fi 


F^ cuts 


AT 


BCDMS (120) 


0.35 - 0.75 


13.25 - 99.00 


59 


32 


1.001 


BCDMS (200) 


0.35 - 0.75 


32.50 - 137.50 


50 


28 


0.998 


BCDMS (280) 


0.35 - 0.75 


43.00 - 230.00 


49 


26 


1.003 


NMC (comb) 


0.35 - 0.50 


7.00 - 65.00 


15 


14 


1.000 


SLAC (comb) 


0.30 - 0.62 


10.00 - 21.40 


59 


59 


0.990 


deuteron 






232 


159 





(b) Number of F2 data points. 



Experiment 


X 


g2, GeV2 


Ff^ 


^2^^ cuts 


J\f 


BCDMS (120) 


0.070 - 0.275 


8.75 - 43.00 


36 


30 


0.983 


BCDMS (200) 


0.070 - 0.275 


17.00 - 75.00 


29 


28 


0.999 


BCDMS (280) 


0.100 - 0.275 


32.50 - 115.50 


27 


26 


0.997 


NMC (comb) 


0.013 - 0.275 


4.50 - 65.00 


88 


53 


1.000 


SLAC (comb) 


0.153 - 0.293 


4.18 - 5.50 


28 


28 


0.994 


nonsinglet 






208 


165 





(c) Number of -^2''^ data points. 



TABLE I: Number of experimental data points (a) F^, (b) F2, and (c) ^2^"^ for the nonsinglet QCD 
analysis with their x and ranges. The name of different data set and range of x and are given in 
the three first columns . The fourth column {F2) contains the number of data points according to the 
cuts: Q2 > 4 GeV2, > 12.5 GeV^, x > 0.3 for i?f and F^ and x < 0.3 for ^3^^. The reduction 
of the number of data points by the additional cuts (see text) are given in the 5th column {F2 cuts). 
The normalization shifts are listed in the last column. 
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-^2,TMc('^'^ ) — / ^ -^2,TMc(-^'Q ) d-X 

Jo 



n + 2 \Q 

+ t ' V i^Kn + 4, Q^) + f ^ ) ,(25) 



2(n + 4)(n + 3) V Q V ^ v ■ ' ^ ■ ^ g2 
where higher powers than (m^/Q^)^ are neghgible for the relevant x < 0.8 region. By 
inserting Eq. fl25|) in Eq. ([5]) we have 

ax ^ 

Ff— '^(x,Q2) = x^(l-x)° ^ e^'^(x) X ^c;")(a,/5)F2^TMc(j + 2,Q') , (26) 

n=0 j=0 

where -^2^^x1^40 + 2,(5^) are the moments determined by Eq. fl25l) . In Fig.([T]) the dashed 
hnes correspond to the NNLO QCD fit adding target mass corrections. 

Despite the kinematic cuts {Q^ > 4 GeV^, = - 1)Q^ + mg, > 12.5 GeV^) 
used for our analysis, we also take into account higher twist corrections to -F|'(x, Q"^) and 
F2{x, Q^) in the kinematic region g2 > 4GeV^4 < < 12.5 GeV^ in order to learn 
whether nonperturbative effects may still contaminate our perturbative analysis. For this 
purpose we extrapolate the QCD fit results obtained for W"^ > 12.5 GeV^ to the region 
> 4GeV^,4 < W"^ < 12.5 GeV^ and form the difference between data and theory, 
applying target mass corrections in addition. Now by considering higher twist correction 
(HT) 

the higher twist coefficient can be extract. Here the operation Otmc[---] denotes taking the 
target mass corrections of the twist-2 contributions to the respective structure function. 
The coefficients h{x, Q"^) are determined in bins of x and and are then averaged over 
Q2. We extrapolate our QCD fits to the region 12.5 GeV^ >W'^> AGeW'^ in Fig.([I]). The 
dashed-dotted lines in this figure correspond to the NNLO QCD fit adding target mass and 
higher twist corrections. There, at higher values of x a clear gap between the data and the 
QCD fit is seen. Figure ([2]) shows the corresponding results for the deuteron data. Figure 
([3]) shows the result of the pure QCD fit for the nonsinglet structure function in NNLO. 
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LO 


NLO 


NNLO 




au 


0.6698 ± 0.0073 


0.7434 ± 0.009 


0.7772 lb 0.009 




K 


3.5104 ± 0.042 


3.8907 ± 0.040 


4.0034 lb 0.033 






0.1990 


0.1620 


0.1000 




du 


1.498 


1.2100 


1.1400 


dv 


ad 


U.boOU it U.Uc)0 


U. i oby it U.U4U 


U. (ooo ± U.U4o 




hd 


3.1685 ± 0.192 


3.5051 lb 0.225 


3.6336 lb 0.244 




Cd 


0.5399 


0.3899 


0.1838 




dd 


-1.4000 


-1.3700 


-1.2152 


AS' MeV 


213. 2± 28 


263.8 lb 30 


239.9 lb 27 


X^/ndf 


538/546 = 0.9853 


523/546 = 0.9578 


506/546 = 0.9267 



TABLE II: Parameter values of the LO, NLO and NNLO nonsinglet QCD fit at = 4 GeV^ The 
values without error have been fixed after a first minimization since the data do not constrain these 
parameters well enough (see text). 



In Table ([IT]) we summarize the LO, NLO, and NNLO fit results without HT contributions 
for the parameters of the parton densities xUy{x,Ql), xdy{x,Ql) and Aq(^. The resulted 
value of x'^/ndf is 0.9853 at LO, 0.9578 at NLO, and 0.9267 at NNLO. Our results for 
covariance matrix for LO, NLO, and NNLO are presented in Table fillip . 



Figure (jlj) illustrates our fit results for xUv{x,Ql), xd^^x 



2\ 



at Qi = 4 GeV^ at NNLO 



with correlated errors. We compare with the results of [30|, 144, |45|] and a very recent 
analysis |.46]. Our results for xUv{x,Ql) and xdv{x,Ql) are in good agreement with the 
other theoretical model at the one a level. 

In Figs. ([5]) and we show the evolution of the valence quark distributions xUy{x,Q'^) 



and xdy{x,Q^) from = lOGeV^ to = 10^ GeV^ in t 



45 



le region x G [10 ^, 1] up to 



46 



471]. With rising values of 



NNLO. We also compared with other QCD analysis [30|, 

the distributions flatten at large values of x and rise at low values. 

Another way to compare the NNLO fit results consists in forming moments of the distri- 
butions Uy{x,Q'^),dy{x,Q'^), and Uj;{x,Q'^) — dy{x,Q'^). In Table HVl we present the lowest 



non-trivial moments of these distributions at 




in NNLO and compare to the respec- 



48 



49| 



tive moments obtained for the parameterizations [30|, 

To perform higher twist QCD analysis of the nonsinglet world data up to NNLO, we 
consider the > 4GeV^,4 < W"^ < 12.5 GeV^ cuts. The number of data points in the 
above range for proton and deuteron is 279 and 278, respectively. The extracted distributions 
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LO 


au 


bu 


ad 


bd 


A Nf=4 
^^QCD 




5.28xl0~^ 










L 

K 


l.boxlU 


1.7oxl0 








ad 


- ( .Oil X iU 


-4.00 X lU 


X.Zo X lU 






Od 


9 fi/1 v 1 n~4 

-Z.D4t X lu 


-L.LL X lU 


D. 10 X lU 


o.O 1 X ±U 




A (4) 


i.yux lu 


Q Q/i 1 n~4 

-o.o4t X iU 


z.oy X lu 


-o. ID X lU 


n T^Q V 1 n~4 


NLO 




bu 


ad 


bd 


A Nf=4 
^^QCD 




8.87x10"^ 










K 


2.39x10 


1.63x10 








0.d 


1 Q/1 V 1 n~4 

-i.o4X iU 


7 QR V 1 n~4 

- ( .oD X iU 


l.Dl X ±U 






Od 


c; 1 n V 1 n~4 

-0. iU X iU 


-4. iy X iu 


s v 1 n~3 

o.oo X iU 


D.U 1 X lU 




A (4) 


o. / i X iU 


r; QQ ^ 1 n~4 

-o.oy X iu 


o.uy X iu 


9 f^7\/ 1 n~4 
Z.O ( X iU 


Q on \/ 1 n~4 

O.OU X lU 


NNLO 


au 


bu 






A Nf=4 
A ' 
^^QCD 


au 


7.61x10^ 












i. 1 O X lU 


1 1 n V 1 n 3 








ad 


-8.41x10"^ 


-6.62x10"'^ 


1.85x10"^ 






bd 


-2.73x10"'' 


-3.73x10-3 


9.79x10-3 


5.98x10-2 




^(4) 


1.08x10"*^ 


-2.74x10"'^ 


1.06x10-^ 


4.19x10-^ 


7.41x10-^ 



TABLE III: Our results for the covariance matrix of the LO, NLO, and NNLO nonsinglet 
= 4 GeV^ by using MINUIT[42| 



QCD fit at 



/ 


N 


NNLO 


BBG 


MRST04 


A02 


A06 




2 


0.3056 ± 0.0023 


0.2986± 0.0029 


0.285 


0.304 


0.2947 




3 


0.0871 ± 0.0009 


0.0871± 0.0011 


0.082 


0.087 


0.0843 




4 


0.0330 ± 0.0004 


0.0333± 0.0005 


0.032 


0.033 


0.0319 


dv 


2 


0.1235 ± 0.0023 


0.1239± 0.0026 


0.115 


0.120 


0.1129 




3 


0.0298 ± 0.0008 


0.0315± 0.0008 


0.028 


0.028 


0.0275 




4 


0.0098 ± 0.0004 


0.0105± 0.0004 


0.009 


0.010 


0.0092 



TABLE IV: Comparison of low order moments from our nonsinglet NNLO QCD analysis at Qq = 4 GeV 



with the NNLO analysis BBG [30], MRST04 



471], A02 



43] and A06 [4 
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for h{x) up to NNLO are depicted in Fig. ([7]) for the nonsinglet case considering scattering 
off the proton target. According to our resuhs the coefficient h{x) grows towards large x. 
Also in this figure HT contributions have the tendency to decrease form LO to NLO, NNLO. 



12| and 



This effect was observed for the first time in the case of fits of A DIS uN data in 

This similar effect was also observed in the fits of F2 charge lepton-nucleon DIS data 



30, |44 



analysis 



, 50, 51|. To compare, we also present the reported results of the early NNLO 
51| in Fig.([7j). Note that the results for h{x) in LO are not presented in the 



44 



BBG model 



30 



5l| . In Ref. 4j], the functional form for h{x) is chosen by 

b 



h{x) 



X 



1 — X 



(28) 



and it is possible to compare h{x) results even in LO. Fig. (El) shows our results for h{x) and 



for the deuteron target up to NNLO. Also we compare the results for the BBG model |51|. 
The same as the proton, HT contributions for the deuteron have the tendency to decrease 
form LO to NLO, NNLO. As seen from Fig. ([7]) and Fig.® h{x) is widely independent of 
the target comparing the results for deeply inelastic scattering off protons and deuterons. 



Our results in low-a; are also in good agreement with 



30 



5l| 



VI. DISCUSSION 



We have performed a QCD analysis of the flavor nonsinglet unpolarized deep-inelastic 
charged lepton-nucleon scattering data to next-to-leading order and derived parameteri- 
zations of valence quark distributions at a starting scale together with the QCD-scale 
Aqcd by using the Jacobi polynomial expansions. 

The analysis was performed using the Jacobi polynomials-method to determine the pa- 
rameters of the problem in a fit to the data. A new aspect in comparison with previous 
analysis is that we determine the parton densities and the QCD scale up to NNLO by using 
the Jacobi polynomial expansion method. The benefit of this approach is the possib ility to 
determine nonsinglet parton distributions analytically and not numerically. In Ref. 52] we 
arrange the MATHEMATICA program to extract Q^) and xd^^x, Q^). 

In this paper the flavor asymmetric combination of light parton distributions x{d — u) 
of Eq. (15) are fixed at = 4 GeV^, as GRS 44| and BBG Q, S applied, and gives 
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a good description of the Drell-Yan dimuon production data 53|]. The first clear evidence 



for the 
CERN 



favor asymmetry of the nucleon sea in nature came from the analysis of NMC at 



54| . In order to have the link with NMC data, we want to study the compatibility 



of the x{d — u) with the NMC result for the Gottfried sum rule (GSR) 55|. This sum rule 
is still actively discussed in problems of deep-inelastic scattering. The GSR, Igsr, can be 
expressed in terms of the parton distribution functions as 

dx 



Igsr{Q') = 



■"2^(3;, Q 



^2 (-^j Q 



X 



1 r 







^ (u,{x, Q2) - d,{x, Q2) j + ^ ( u{x, Q2) - d{x, Q2) 



dx 



1 2 
3 ^3./o 



u{x, Q^) — d{x, Q"^) \ dx . 



(29) 



In the derivation of the above equation, the asymmetry of nucleon sea was assumed. The 
NMC measurement |15^ implies at = 4 GeV^ 



(c/(x, Q^) - u{x, Q^))dx = 0.148 ± 0.039 



(30) 



which was the first indication that there are more down antiquarks in the proton tha n up 



27| 



antiquarks. On the other hand this value is reported 0.118 ± 0.012 at = 54 GeV^ 
Now it is interesting to obtain this value for the parametrization of Eq. (|T5l) which we used 
in our QCD analysis. By integration of this distribution we obtain ^ 0.1 which is smaller 
than the reported results in the literature. However, the NMC Collaboration gives the Igsr 
experimental value at = 4 GeV^ 54 1 



J^^j'^(g2 = 4 Qgv2) = 0.235 ± 0.026 . 



(31) 



By using Eq. ( l29|l we obtain the GSR value about 0.267 with which the existing measure- 
ments are almost compatible within error. It seems that although the value of jQ{d — u)dx 
is smaller than the values in the literature, the parametrization of Eq. (fT5|) can give a good 
description of the E866 experimental data '2^- Also we should notice that the GSR does 
not belong to the strict sum rules in QCD and it is necessary to receive not only QCD 



58 



59 



60 



6l| 



corrections but anomalous dimensions as well [S^l, [s?. 

Now it is interesting to compare the NNLO theoretical QCD theoretical prediction for 
;he Gottfried sum rule |58|1 with NMC data. The recent step in this direction was done in 



62] • According to this paper we add the QCD two-loop correction to the Gottfried sum 
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rule and we refine the GSR value to about 0.12%. Also we obtain the value of /G5ij(0.004 < 
X < 0.8,4 GeV^) = 0.267 which is well compatible with the neural parametrization results, 
e.g. 0.2281 ± 0.0437 ^ within errors. 

In the QCD analysis we parameterized the strong coupling constant in terms of four 
massless fiavors determining Aqcd- The LO, NLO, and NNLO results fitting the data, are 

AgJI^s = 213.2 ± 28 MeV, LO, 

A'^^^^ = 263.8 ± 30 MeV, NLO, 

Agg^s _ 239.9 ±27 MeV, NNLO, (32) 

These results can be expressed in terms of q;s(M|): 

a,(M|) = 0.1281 ±0.0028, LO, 
a,(M|) = 0.1149 ± 0.0021, NLO, 

a^(M|) = 0.1131 ±0.0019, NNLO. (33) 

Note that in above results we use the matching between Uf and fiavor couplings cal- 



63| . To be capable to compare with other measurement of Aqcd we adopt 



culated in Ref. 
this prescription. 

The as(M|) values can be compared with results from other QCD analysis of inclusive 
deep-inelastic scattering data in NLO 



A02 


[48] 




=0.1171 ±0.0015 


ZEUS 


[64] 




=0.1166 ±0.0049 


HI 


m 




=0.1150 ±0.0017 


BCDMS 


[36J 




=0.110 ±0.006 


GRS 


m 


«s(M|) 


=0.112 


CTEQ6 


[65J 


«s(M|) 


=0.1165 ±0.0065 


MRST03 


[66] 




=0.1165 ±0.0020 


BBG 


m 


a.(M|) 


=0.1148 ± 0.0019 


KK05 [67] 


a.(M|) 


=0.1153 ± 0.0013(stat) 



±0.0022(s?/st) ±0.0012(norm) 



BE (pol)[68j: 
AK (pol) [4]: 



a,(M|)=0.113 ±0.004 
a,(M|)=0.1141 ± 0.0036 
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The NNLO values of as(M|) can also be compared with results from other QCD analysis 



A02 [48] 
GRS \U\ 
MRST03 [66] 
SYOl(ep) [69] 
SY01(z/N) [69] 
A06 [49.1 



a,(M|)=0.1143 ±0.0014 
«,(M|)=0.111 
a,(M|)= 0.1153 ±0.0020 
a,(M|)=0.1166 ±0.0013 
a,(M|)=0.1153 ±0.0063 
a^(M|)=0.1128 ± 0.0015 



BBG [301: 



a,(M|)=0.1134 



BM07 [70]: 
KPSOO(i/iV) [JJ: 



KPS03(i/A^) [15 



+0.0019 
-0.0021 

a,(M|)=0.1189 ±0.0019 

a,(M|)=0.118 ± 0.002 {stat)± 0.005 (syst) 

± 0.003 (theory) 

a,(M|)=0.119 ± 0.002 {stat)± 0.005 (syst) 

± 0.002 (threshold) tofot (scale) 



and with the value of the current world average 



0.1189 ±0.0010 



which has been extracted in 



7l| recently. 



(34) 



We hope our results of QCD analysis of structure functions in terms of Jacobi polynomials 
could be able to describe more complicated hadron structure functions. We also hope to 
be able to consider the N^LO corrections and massive quark contributions by using the 
structure function expansion in terms of the Jacobi polynomials. 
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FIG. 1: The structure function i^f as a function of in intervals of x. Shown are the pure QCD fit 
in NNLO (solid line) and the contributions from target mass corrections (dashed line) and higher twist 
(dashed-dotted line). The vertical dashed line indicate the regions with > 12.5 GeV^. 
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FIG. 2: The structure function as a function of in intervals of x. Shown are the pure QCD fit 
in NNLO (solid line) and the contributions from target mass corrections (dashed line) and higher twist 
(dashed-dotted line). The vertical dashed lines indicate the regions with W'^ > 12.5 GeV^. 
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FIG. 3: The structure function F^^^ as a function of in intervals of x. Shown is the pure QCD fit 
in NNLO (solid lines). 
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FIG. 4: The parton densities xu^ and xd^ at the input scale Qq = 4.0 GeV^ (solid line) compared 



451], 



with results obtained from NNLO analysis by BBG (dashed- line) A05 (dashed-dotted line) 
MRST (dashed-dotted-dotted line) [46], and GRS (dashed-dashed-dotted line) j44|. The shaded areas 
represent the fully correlated one a statistical error bands. 
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FIG. 7: The higher twist coefficient h{x) for the proton data as a function of x up to NNLO (solid 
ine) compared with results obtained by GRS (dashed line) and BBG (dashed-dotted line) 
51|. 
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FIG. 8: The higher twist coefficient h{x) for the deuteron data as a function of x up to NNLO 



(soHd line) compared with results obtained by BBG (dashed-dotted line) 
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